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On the Velocity of Steady Fall of Spherical Particles through 

Fluid Medium, 

By E. Cunningham, M.A., Fellow of St. John's College, Cambridge, 
Assistant Professor of Applied Mathematics, University College, London. 

(Communicated by Sir Joseph Larmor, Sec. K.S. Eeceived January 18, — Eead 

January 27, 1910.) 

1. Eecent investigations^ have called attention to the question as to how 
closely the rate of fall of particles through a fluid when very small approxi- 
mates to that given by the law obtained by Stokes from hydrodynamical 
reasoning. This paper deals theoretically with two of the main sources of 
divergence from that law. 

Firstly, taking the particles to be smooth spheres moving through a gas, it 
is shown that the deviation to be expected on account of the diameter of the 
particle being small, compared with the molecular free path, is extremely 
small for such particles as have been experimented with ; and a modified 
formula is given which may be taken to hold approximately for lower 
pressures or particles of smaller dimensions. Secondly, an approximate 
treatment is made of the effect of the simultaneous presence of a large 
immber of particles moving with the same velocity through the fluid, and it 
is found that the force required to maintain the motion of one of them 
depends not only on the diameter, but on the ratio of the diameter to the 
distance between the particles, and in such a way as to increase rapidly as 
this ratio increases beyond the value 01. 

With the figures obtained by J. J. Thomson,f it is shown that these 
corrections together would imply an increase of not more than 0*07 x 10"^'^ 
in the charge on the gaseous ion. 

'^' See the preliminary experiments of Zeleny with spores (British Association, 1909)^ 
and of Buller ('Nature,' April, 1909); also the later concordant results of Zeleny with 
spherical particles ('Nature,' Dec, 1909); also Millikan ('Physical Review,' Dec, 1909). 
Since this paper was prepared, Prof. Larmor has informed the writer that Prof. Zeleny is 
pursuing the subject of fall in gases at diminished pressures, with a view to throwing light 
on the nature of a molecular collision with a solid body, such as is required for the theory 
of the radiometer and other purposes {of. Maxwell, ' Phil. Trans.,' 1879). For convenience 
of reference, a short calculation on this point has been added. The formula obtained in 
the first part of this paper agrees with the general conclusion that experiments of this 
kind cannot lead to a determination of molecular magnitudes, as at first sight seemed not 
impossible. 

t J. J. Thomson, 'Phil. Mag.,' vol. 5, 1903. 
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On the Bate of Diffusion of a Cloncl of Particles through a Gas. 

2. Suppose that a cloud of particles of mass M at a concentration of 
N" particles per unit volume is diffusing through a gas, the mean velocity of 
the particles being Vo, the mass of the gas molecule being m, and the number 
of molecules per unit volume being n. It is desired to find the mean force 
per particle required to maintain tlie velocity Vq. 

It will be assumed in the first instance that the collisions are of the nature 
of impulses between smooth elastic spheres, and that a collision takes place 
between a molecule and a particle when the distance apart of their centres 
is a. 

Consider the circumstances of collisions of this type with a single particle, 
the centre of the molecule at collision being within a small area 

dco = a^ sin d9 d^ 

of the sphere of radius a surrounding the particle, (r, 0, (p) being spherical 
polar co-ordinates having the direction of Vo as axis. 

Let the velocity of the molecule be (v., v, w) and that of the particle 
(Vo + U,V,W). . 

The probability of molecular velocity lying within a range dudvdw 
about (tt, Vy to) is 

where A (Tr/hm)^ = 1 ; 

and that of the velocity of the particle lying within a range dJJ dV dW is 

Be~^^ (ir^+ v^+ w^-) ^jj ^y ^^^ 

where B {irpiKf = 1. 

In order that a collision may take place within the given limits in a given 
interval of time ht, the centre of molecule must at the beginning of it lie 
within a cylindrical volume of base dco and height (U + Vocos^— 2/^) ht, with 
the restriction that (IT + Vo cos — ^^) must be positive. 

The magnitude of the impulse I is (2Mm/M4-m) (U + Yocos0 — 2^). 

Thus the probable impulse in the direction of Vo on the particle per unit 
time is 



f 




^^._/,(^,2+^.2+,,2) J ^^^ Q ^^^ ^jj _^ y^ ^^g 0—v) d<o d.u dv dw, 



the limits for v, w, 0, (p being unrestricted, but not so for m. This expression 
is equal to 



. TT o 9 2M??^ 
hrrb M + 7n 



L-hmu^ (U + Vo cos e-^uf sin cos du d0. 
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Multiplying by the probability of the given velocity (U, V, W) and 
integrating, the mean impulse per particle per unit time becomes 

I AB^i(|Y47r|^j^'| f f fc-M-«^^+Mu«) (tn- Vo cos O-uf sin 6 cos dO du d\], 

the limits for 6 being and tt, II and % taking all values for which 
'U < U + Vo cos 6, 

The evaluation of the last integral is only possible on the assumption that 
Vo is small compared with the mean relative velocity. Making this 
assumption, and calling the integral K, we have 

-where a =U-~i^, ^ = (MU + m'ie)/(M + m). 

Thus 



.^ ^co 



K= \ rTTFT— \r 1 I e-^""^'^'/^^+^^>(a + VoCos^)2sin6>cos6>^6>6?a. 
If the lower limit for a were 0, the value of this expression would be 
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TT I 2 



h (M + m) 







'The error in taking this value for K is 



TT 



1 



.7?.(M + m). 
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-and is therefore of the order A^o/^ times K, where c is the mean value of a, so 
that to the first order the above value may be taken for K. 
Hence the mean impulse per particle per unit time is 

{ABn {irjhf 47ra-V(M + m)} f Vo v" { ^r (M + m)imn'M.^} 
= f a•2Vo^lv'{7^M7?^/(M^-m)2}.* 

This, then, is the force per unit particle required to maintain a drift Vq of 
the particles relative to the gas. It represents the rate of transfer of 
imomentum from a particle to the surrounding gas by reason of collisions. 

It may be expected that, in the case where the particles are few in number 
•compared with, and of the same order of magnitude as, the molecules, the mass 
velocity of the gas in the neighbourhood of a particle will be very slightly 
affected by the motion of the particle ; but if the diameter of the particle 
becomes comparable with the mean free path of the molecule, it will become 

* This expression is also obtained as a particular case of a general formula by Langevin 
■('Annales de Chimie et de Physique,' vol; 5, 1905, p. 266), and is applied by him to the 
consideration of the mobility of ions. 
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impossible to neglect the mass velocity produced in the gas by the motion of 
the particle. It will be necessary to substitute for Vo in the above expression 
for the transfer of momentum^ the difference between the mean velocity of 
the particle and the mass velocity of the gas. 

It may be remarked that though a cloud of particles has been spoken of 
above, the result will hold equally well for the mean velocity of a single' 
particle taken over a time sufficiently long to include a very large number of 
collisions. This is the case that must be considered in examining the 
deviation from Stokes' Law, which is demonstrated only for a single particle- 
in an infinite fluid. 

3. We may obtain now a modification of the law, taking into account the 
machinery by which the transfer of momentum from particle to fluid takes 
place. This is not considered in the ordinary hydrodynamical investigation « 
Instead, a Icinematic assumption is made, namely, that the mass velocity in. 
the neighbourhood of the particle is Y. Let us assume, in place of this, that 
the mass velocity is hY, Then the rate of dissipation of momentum on the 
hydrodynamical theory is QTr/juakY, a being used now for the effective radius 
of the particle, inasmuch as it is assumed to be large compared with the 
molecule. In the steady motion this diffusion of momentum must just 
balance the momentum communicated by the particle to the gas. Taking 
M large compared with m, this is, from the above formula, 

Hence 6 'Kixak = f a / (-7- j ^^^^ (V— Z;V), 



or 



7 _ 4 an (irmjlvf — /1 _l '^ ^ /3 tt' 



4 an {Trmjli)^ -\- pTr/jb \ 4: a v 2 

where P is the mean free path of the molecules in the air. Hence the force 
required to maintain the velocity V is 

X = 6 TTfiaJcV = 6 TTfiaY (1 + 1*63 l/a)~^ ; 

that is, the limiting velocity of a falling particle will exceed that given by 
the ordinary formula in the ratio of (1-f l*63Z/a) to unity. This ratio for a 
given gas should therefore differ from unity by an amount which varies 
inversely as the density. 

For air at ordinary pressure and temperature / may be taken to be about 
0*9 X 10"^, so that the limiting velocity would, on the present calculation, be 
double the usual value for a particle of radius about 1*5 x 10"^ For particles 

■^ I is the mean free path obtained from the viscosity by the formula /x — J pcl^ not the 
value corrected for persistence of velocity {vide Jeans, ' Kinetic Theory of Gases,' p. 250). 
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•of the magnitude of the smallest spores in Zeleny's earlier experiments the 
ratio would be about 1*07. 

The nature of the collision between molecule and particle has, so far, been 
assumed to be that which occurs between smooth elastic spheres. If this is 
replaced by supposing that each molecule became entangled in the particle, 
so that on re-emission there was no correlation between its velocity and that 
with which it collided, a similar calculation can be easily made. The impulses 
by which the molecules are thrown out from the particle have no resultant, if 
all directions of re-emission are supposed equally probable. The impulse at 
•collision in the direction of the velocity Yi) is 

.^^1 {Vo-f(U-'^0 cos 6>-f(V-^/;) sin 6}. 

If the integration be effected as before, the mean impulse is found to bear to 

(1 -f cos^ ^) sin ^ ^/0 to 2co^^ 6 ^\n 6 d6, 

Jo 

that is, it is exactly double. 

If, now^ it is assumed, following Maxwell, that a fraction "/" of the colliding 

molecules experience collisions of the earlier type, and the remainder of the 

latter type, the equation giving the limiting velocity becomes 

Maxwell {loc. cit.) deduces from experiments of Kundt and Warburg that 
for impacts of air molecules on glass the fraction / is approximately half. 
This would reduce the divergence from Stokes' Law to two-thirds of that 
obtained by putting / = 1. If no molecules experience perfect reflexion, 
/ would be zero, and the divergence would be reduced to one-half. If it be 
assumed that for the majority of material particles / lies between the values 
■0 and ^, the result is that Stokes' formula will yield a velocity equal to half 
the actual for particles whose radii lie between 0*75 x 10"^ and 1 x 10"^, while 
for particles of diameter 2 x 10"*, such as those used in Zeleny's earlier 
experiments, it will yield a velocity between 1*05 and 1*036 too small. 
Experiments at lower pressures should show greater deviations, and may 
throw light on the value of /. 

On the M'lfMtal Injlmence of the Fartictes in a Cloud. 
4. A possible source of deviation from Stokes' Law for a cloud of particles 
of dimensions for which the hydrodynamical condition of no surface slipping 
is satisfied, is the mutual influence of the particles, and the effect of the 
boundaries. The exact solution of this problem is clearly beyond analysis. 
But an approximate method of attack suggests itself, based on the fact 
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that if a cloud of particles is uniformly distributed throughout aii- 
enclosed volume of fluid and are moving with equal parallel velocities, the- 
drift of the fluid with each particle in the direction of its motion in the 
immediate vicinity of the particle must be compensated by a drift in the 
opposite direction between the particles. In fact, if the distribution is 
uniform the mean flux of fluid across any area is zero if the fluid is 
contained within fixed boundaries. In Stokes' problem this condition 
is absent. 

The method employed to obtain an idea of the magnitude of the 
influence of this consideration is to imagine that each particle is effectively 
limited to moving within a concentric spherical mass of fluid whose boundary 
is fixed. If, for instance, supposing the particles to be of radius a, and 
surrounded by concentric spheres of radius 5, the arrangement of particles 
being such that these latter spheres are packed in the closest possible 
arrangement, the fluid in the spaces between them being supposed occupied 
by fluid at rest, the conditions of the problem are to some extent realised ; 
but the velocity of each particle in its own sphere will be somewhat less 
than the actual limiting velocity owing to the restriction imposed on the 
fluid between the spheres. We may therefore expect by this method to- 
obtain a lower limit for the velocity sought in terms of the radius of each 
particle and the number of particles per unit volume. 

The problem is therefore to be solved of the drag due to viscosity on 
a particle moving with velocity V at the centre of a s|)herical envelope. 
This can be effected exactly, the method adopted being that given by Lamb* 
for the solution of Stokes' problem. 

Taking the axis of x in the direction of motion, the velocity at a point can 
be expressed in terms of spherical harmonics of the first order. 

Putting pi — ocw, p~2 = Axjr'^, 



^1 = j3x, (f)-2 ~ Bx/ 



,3 



and using the formula of type 



u 



__ 1 ^ r___r^____ Sgw, fir/2n + 3 ^ _Pn_\ , V ^<l> 



fi L2 (2 n + 1) dx (n + 1) (2 7i + 1) (2 n + 3) dx r^'^+i J ax 

we obtain 

^ r r2-3a;'n ^2A. fl r^-^S^n , j,r^-3x^ , ^ 

^■^'— T~ i ^ + r r + Q— i J-T~^ i" + ^ ?— +^^ 

o/x L 5 J 6 pi [r 4r'^ j r^ 



a 



, XX. (Ay I J ty -w-y X (f 

^ — — — — xy + . —6 ij-~ ^ 

OL Jlx. Jb/i> fy -ry XZ 

W = — — - — ' XZ -T -r: 5 — 0x5 — ^ . 

10 fjb 2 fjb r^ T^ 



* 'Hydrodynamics,' 3rd edit., § 324-5. 
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!N"ow at the surface of the particle where r — a, 

u = Y, V = 0, w = 0, 

and at the sphere r = b u z= v = lo — 0, 

These conditions give the equations 

o?a ,2A,o ^r a, A 3B r. 

K) fjb S fia 10 /x 2 fjbcr fjp 

The resultant traction on the sphere is found to be, exactly as in the 

simpler problem,* 

X= -47rA. 

The value of A obtained from the four equations above gives 



a result which reduces at once to Stokes' law when hja is large. 
Expanded in powers of (a/&), the result is 

It was said above that the value so found is a lower limit for the actual 
velocity under a given force, owing to the conditions imposed being possibly 
more rigid than the problem demanded. 

5. If now round each particle be described a sphere of radius h\ such that 
the new spheres overlap, but only by so much that they just include the 
whole of the fluid, each particle being again supposed to move in its own 
sphere, a new value for the velocity will be obtained. The arrangement of 
particles being as before, V will be the distance from a corner to the centroid 
of a tetrahedron whose edge is 2 &, so that V =z h ^ (3/2). 

It is not clear, however, that in this way an upper limit to the velocity of 
the particles is obtained, though the following considerations point in that 
direction. 

In either case the argument is that if the velocity system for the fluid be 
obtained when the fluid within one of the spheres is set in motion by the 
enclosed particle, the external fluid being at rest ; and if the velocity 
systems thus obtained for all the particles are superposed, the resultant 
velocity system will approximate to the actual in the sense that in the 
neighbourhood of each particle the conditions imposed by the continuity of 

^ Loe. cit.^ p. 552, equations (5) and (13). 
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velocity is satisfied, and that the hydrodynamical equations, being linear, are 
also satisfied. 

Thus, in the former case, the region external to all the spheres of radius h 
is kept at rest ; while, in the latter, the motion in a region internal to two or 
more spheres is compounded of the contributions from the motions in the 
separate spheres. 

But the difficulty is thus introduced, that, at the surfaces of the auxiliary 
spheres, discontinuity occurs in the velocity gradient. In the actual problem 
the discontinuity is smoothed out. 

The diagram below shows, for hja = 5, the variation of the vertical velocity 
of the fluid as obtained above : (i) along a line through the centres of two of 




Ai Aa represents velocity of particles. 

base line for curve (i). 



the spheres supposed at distance 2 h apart in a horizontal line ; (ii) along a 
median of an equilateral triangle of side 2 5, the^-e being supposed three 
spheres with their centres at the corners. The continuous line refers to the 
velocity found by the use of the larger spheres, the broken line to that given 
by the smaller. Comparison of these curves and of others not given here 
indicates fairly clearly that the former is more likely to be near to the true 
facts. Accordingly the values are tabulated below of the ratio of the forces 
required to maintain a velocity V according to the method adopted above, 
using the larger sphere, and, according to Stokes' Law, for different values of 
h' Ja ; the second row of the table gives the corresponding value of hja. 
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6. Assuming, now, that the modified laws obtained above are correct^ we 
may see what correction . they imply as necessary to J. J. Thomson's calcula- 
tion of the charge on the gaseous ion. 

The figures given by him fall into four groups corresponding to values of 
the charge on the ion— 

101^6 = 3-2, 3-50, 3-45, 3-39, 

respectively. Taking the figures given for the number of drops per cubic 
centimetre, we are able to calculate the corresponding values of ¥ Ja, and 
thence the correction to be made as in § 5. The result is that e must be 
increased by 6*2, 7*4, 6, 5*9 per cent, in the four cases respectively. For the 
calculated values of the radii of the drops, however, the correction on account 
of the first cause (§ 3) amounts to about 4 per cent, and is in the opposite 
sense, and therefore corresponds to a diminution of 4 per cent, in e. 

On account of both causes we may say that e is to be increased by about. 
2 per cent., that is by about 0*07 x lO""^^. 

It is clear, therefore, that neither of the causes discussed accounts for more 
than a small fraction of the discrepancy between the older and the more- 
recent estimates.* 



* Vide Rutherford, Presidential Address to British Association, Section A, 1909 ; alsov 
Millikan, ' Physical Eeview,' Dec, 1909. 
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